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We consider transport in a three terminal device attached to one superconducting and two nor- 
mal metal terminals, using the circuit theory of mesoscopic superconductivity. We compute the 
nonlocal conductance of the current out of the first normal metal terminal in response to a bias 
voltage between the second normal metal terminal and the superconducting terminal. The nonlocal 
conductance is given by competing contributions from crossed Andreev reflection and electron co- 
tunneling, and we determine the contribution from each process. The nonlocal conductance vanishes 
when there is no resistance between the superconducting terminal and the device, in agreement with 
previous theoretical work. Electron cotunneling dominates when there is a finite resistance between 
the device and the superconducting reservoir. Dephasing is taken into account, and the charac- 
teristic timescale is the particle dwell time. This gives rise to an effective Thouless energy. Both 
the conductance due to crossed Andreev reflection and electron cotunneling depend strongly on 
the Thouless energy. We suggest to experimentally determine independently the conductance due 
to crossed Andreev reflection and electron cotunneling in measurement of both energy and charge 
flow into one normal metal terminal in response to a bias voltage between the other normal metal 
terminal and the superconductor. 



I. INTRODUCTION 

Crossed Andreev reflection^ transforms an incident 
electron from one conductor, attached to a supercon- 
ductor, into a hole in a geometrically separated second 
attached conductor. In an alternative, equivalent pic- 
ture, two quasiparticles from two separate conductors 
are transferred into a superconductor as a Cooper pair. 
Electrons can also be transferred between the conduc- 
tors by electron cotunneling, where an incident electron 
tunnels via a virtual state in the superconductor. The 
nonlocal conductance, defined in a three terminal device 
as the current response in one normal metal terminal 
to a voltage bias between the other normal metal and 
the superconductor, is determined by contributions from 
both crossed Andreev reflection and electron cotunnel- 
ing. Crossed Andreev reflection and electron cotunneling 
give opposite contributions to the nonlocal conductance. 
In this way, crossed Andreev reflection competes with 
electron cotunneling. The realization of a system where 
crossed Andreev reflection can be observed, has been the 
aim of both experimental^ and theoretical 5 -! 6 -^^ work 
lately. This interest is due to the fact that crossed An- 
dreev reflection is an inherently mesoscopic phenomenon, 
with the prospect of creating entangled electrons. 10 ' 11 ! 12 

The nonlocal conductance of a device where two nor- 
mal conductors are tunnel coupled to a bulk supercon- 
ductor was calculated in second order perturbation the- 
ory for quantum tunneling^ in Ref. [H The conductance 
originating from crossed Andreev reflection was predicted 
to exactly cancel the conductance due to electron cotun- 
neling. Subsequently disorder— and higher order quan- 
tum interference effect s 15 ' 16 have been incorporated into 
this approach, and the noise and cross correlations have 
been considered^ Ferromagnetic contacts were also con- 



sidered in Ref. [f| Crossed Andreev reflection is favored in 
an antiparallcl configuration of the magnetizations, since 
Cooper pairs in singlet superconductors consist of two 
electrons with opposite spin. Electron cotunneling is fa- 
vored in a parallel configuration. 

The predicted dependence of the nonlocal conductance 
on the magnetization configuration was observed ex- 
perimentally in a hybrid superconducting/ferromagnetic 
device^ Subsequently, a bias dependent nonlocal conduc- 
tance was observed in a more complicated geometry with 
only normal metal contacts to the superconductor^ For 
bias voltages corresponding to energies below the Thou- 
less energy associated with the distance between the two 
normal terminals, a nonlocal signal with sign correspond- 
ing to electron cotunneling was seen. Thus, in contrast 
to the results of Refs. [H and [lj, experiments showed a 
finite nonlocal conductance at low bias. Additionally, the 
sign of the nonlocal signal in Ref. changes when the 
bias voltage exceeds the Thouless energy, and this was 
interpreted as a consequence of crossed Andreev reflec- 
tion dominating the nonlocal signal. These experimental 
findings arc currently not understood. 

In previous theoretical works, it is assumed that su- 
perconducting properties, e.g. the magnitude of the gap, 
are not modified by the presence of the contacts. This 
assumption is valid as long as the coupling between the 
normal/ferromagnetic conductor and the superconductor 
is weak or has a small cross section compared to the su- 
perconducting coherence length. None of the mentioned 
theoretical works describe a dependence in the conduc- 
tances on the Thouless energy. 

The circuit theory of mesoscopic transpor t 18 ' 19 is a 
suitable framework to understand transport properties 
of mesoscopic small normal metal/superconducting sys- 
tems. Circuit theory is a discretization of the quasiclassi- 
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cal theory of superconductivity^ and can treat nonequi- 
librium effects and dephasing. A circuit is modeled as a 
network of cavities, connectors, and terminals - similar 
to the way we understand classical electric circuits. Ter- 
minals are voltage probes in local thermodynamic equi- 
librium, whereas cavities can be driven out of equilib- 
rium. Cavities and terminals may be normal metals or 
superconductors. The connectors can represent physical 
interfaces between cavities and terminals or model dif- 
fusion. Connectors representing interfaces are described 
by their sets of transmission probabilities. "KirchhofFs 
rules" determine the matrix Green's functions (poten- 
tials) of the cavities and the matrix currents through the 
connectors. The matrix currents describe not only the 
flow of charge and energy, but e.g. also the flow of quasi- 
particlc correlation. Circuit theory has been applied suc- 
cessfully to explain various phenomena in superconductor 
and normal metal/ferromagnet hybrid structures, like the 
proximity effect ji 8 . multiple Andreev reflections j21 spin 
transport) 22 ' 23 and unconventional superconductivity^ 
Circuit theory has also proved to be a successful ap- 
proach to calculate the full counting statistics of hybrid 
structures. 25 ' 26 ' 27 ' 28 A circuit theory of magnetoelectron- 
ics has been developed for hybrid systems of fcrromagncts 
and normal metals, see Ref. and references therein. 

We apply circuit theory to calculate the nonlocal con- 
ductance of a three terminal device with contacts to one 
superconducting and two normal metal terminals. We 
substantially generalize previous theoretical approaches 
by computing analytically the nonlocal conductance with 
general connectors ranging from e.g. ballistic point con- 
tacts via diffusive contacts to tunnel junctions. We take 
the proximity effect into account, in which superconduct- 
ing correlations affect the spectral properties of a normal 
metal. We also take dephasing into account, where the 
dwell time gives rise to an effective Thouless energy. We 
consider low bias transport so that the relevant energies 
are much smaller than the superconducting gap. The 
model we consider has a simple and generic geometrical 
structure. We do not consider Josephson effects. We re- 
cover several aspects seen in experiments^ Crossed An- 
dreev reflection and electron cotunncling do not cancel 
each other. However, in the limiting case of strong cou- 
pling between the device and the superconducting ter- 
minal, our results agree with previous theoretical work, 
and the nonlocal conductance vanishes. The differen- 
tial conductances depend on the effective Thouless en- 
ergy. A dependence on the Thouless energy has been 
observed experimentally^ However, we do not find an 
agreement with the sign of the measured nonlocal con- 
ductance above the Thouless energy for the simple model 
we study. 

The paper is organized in the following way: In Sec- 
tion HU we give an overview of our model and the circuit 
components. In Section HTT1 we present the mathematical 
description and calculation that determines the conduc- 
tances associated with the various transport processes in 
the system. We show numerical results for some experi- 



mentally relevant systems in Section HVl Finally, we give 
our conclusions in Section fVl 



II. MODEL 




FIG. 1: Our circuit theory model. A normal metal chaotic 
cavity (c) is connected to one superconducting (S) and two 
normal metal terminals (Ni and N2). The three connectors 
are described by their sets of transmission probabilities. A 
coupling to ground represents the "leakage of coherence" (see 
text). 

We consider a three terminal system where one super- 
conducting terminal and two normal metal terminals are 
connected to a small normal metal cavity. We assume 
that the cavity is large enough that charging effects can 
be neglected, and that the Green's functions are isotropic 
due to chaotic scattering. A physical realization of the 
chaotic cavity could be a small piece of diffusive metal, 
embedded in a circuit by e.g. tunneling contacts to the 
terminals. The assumptions on the chaotic cavity are 
quite general and can also be fulfilled for e.g. a quan- 
tum dot with ballistic point contacts if the conductance 
of the contacts is much less than the Sharvin estimation 
of the cavity conductance^ The circuit theory repre- 
sentation of our model is shown in Fig. [T] The normal 
terminals Ni and N2, and the superconducting terminal 
S, are connected to the chaotic cavity c through general 
connectors represented by their sets of transmission prob- 
abilities where i = 1, 2, S, and the index n num- 
bers the conductance channels. These connectors can 
represent anything from ballistic point contacts to tun- 
nel junctions.— For a ballistic connector all transmission 
eigenvalues are equal to 1 for the propagating channels 
and otherwise. For a tunnel junction, all transmis- 
sion probabilities are small. Dephasing is represented in 
the circuit diagram by a coupling to ground, although 
no energy or charge current can flow to this terminal. 
The dephasing will be discussed in more detail in Sec- 
tion Uni Our model has a generic geometrical structure 
and will capture the physics of crossed Andreev reflection 
and electron cotunncling for a wide range of systems. 
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Let us now identify the various transport processes in 
the system. We expect the following contributions to the 
current: Electron cotunncling (EC) between terminals 
Ni and N2, direct Andreev reflection (DA) between the 
superconductor and either normal terminal Ni or N2 , and 
crossed Andreev reflection (CA) between the supercon- 
ductor and both normal metal terminals Ni and N2. In 
direct Andreev reflection, an injected particle from one 



terminal gives rise to a reflected hole in the same ter- 
minal, whereas in crossed Andreev reflection an injected 
particle from terminal N2 (Ni) gives rise to a reflected 
hole in terminal Ni (N2). These processes are illustrated 
in Fig. 

Semiclassical probability arguments show that the 
spectral charge current in the connector between Ni and 
c at energy E has the following structure^ 
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(a)Crossed Andreev reflection: A particle from N2 with 
energy eV-2 and a particle from Ni with energy —eV^ form 
a Cooper pair in S. 



eV 2 
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(b)Electron cotunneling: A particle from N2 at energy eV2 
tunnels through the cavity c into Ni . The density of states 
in the cavity c is suppressed due to the proximity effect 
from the superconducting terminal. 



FIG. 2: Transport processes in the three terminal device. 



h{E ) = {f 2 (E) - h(E)] + [1 - h{E) - h{-E)\ + {1 _ h{E) _ h{ _ E)] t (1) 



where /j (±E) denotes the Fermi-Dirac distribution func- 
tions in normal terminal i at energy ±E. The energy 
dependent conductances G(E) are even functions of en- 
ergy. Andreev reflection couples an electron with en- 
ergy E in terminal Ni to an electron with energy —E 
in either terminal Ni (DA) or N 2 (CA). The factor 
2 for direct Andreev reflection takes into account that 
two charges are transmitted in this process. We di- 
vide currents and distribution functions into even and 
odd parts with respect to energy The even part con- 
tributes to spectral charge current, and the odd part 
contributes to spectral energy current. We therefore con- 
struct I T ,i(E) = {h(E) + h(-E)\ /2: 



It,i(E) — - ( -Gec + tGca 
e \ 2 2 



2Gda Iit,i 
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where /ix,i are the reservoir distribution functions deter- 
mined by Fermi-Dirac functions^ The total charge cur- 
rent is found by integrating over all energies /charge, 1 = 
/ AEI'y i. We define the nonlocal differential conduc- 
tance as the current response in one normal terminal to 
a voltage between the other normal terminal and the su- 



perconductor. Using Eq. this quantity becomes 



dl, 



charge,! _ 



-JdE[G EC (E)-G CA (E)] dfiE dE eV2 \ 

(3) 



where V2 is the voltage in terminal N2. At low tem- 
perature, the derivative of the Fermi function gives a 
delta-function at energy eV^. The integral then gives 
d/charge,i/0V2 = G-®c{pV'i) ~ G C A(eV 2 ), thus the nonlo- 
cal differential conductance is determined by the differ- 
ence of the cotunneling conductance and the crossed An- 
dreev reflection conductance for quasiparticles at energy 
eV%. Consequently, a measurement of the nonlocal con- 
ductance does not uniquely determine the conductance 
for both processes. 

In Section lTnl we show that according to the circuit the- 
ory, Gec — Gca is positive at all energies for the generic 
network considered. This means that the conductance 
resulting from electron cotunneling is larger than the con- 
ductance resulting from crossed Andreev reflection, and 
thus the nonlocal differential conductance is always pos- 
itive. 

The odd part of the current in Eq. ([T]) with respect 
to energy contributes to energy transport. Direct An- 
dreev reflection does not contribute since two particles 
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of energies E and —E arc transmitted. The spectral en- 
ergy current in the connector between Ni and c has the 
following structure 40 



i and the cavity i: 
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The total energy current is obtained from the spectral 
energy current 7 L)1 by / C ncr g y,i = / dEEI hA /e. The 
energy current into the terminal is related by the heat 
capacity to the rate of change of the temperature. Thus 
a nonlocal differential conductance for energy transport 
could in principle be measured by considering the heat 
flow into terminal Ni . We define the nonlocal differential 
conductance for energy transport as the energy current 
response in one terminal to a voltage between the other 
normal terminal and the superconductor. From Eq. (j4|) 
this quantity becomes 
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dV 2 



dE-[G EC (E) 
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Gca(E)} 



df(E-eV 2 ) 



dE 



(5) 



At low temperatures, this gives dI enelgyt i/dV2 = 
V2[GEc(e^) + GcA(eV2)], thus the nonlocal differential 
conductance for energy transport is determined by the 
sum of the cotunncling conductance and the crossed An- 
dreev reflection conductance for quasiparticlcs at energy 
eV 2 . 

The discussion above shows that the conductance 
of electron cotunncling and crossed Andreev reflection 
can be determined independently from two experimen- 
tal quantities. Measurements of the nonlocal differential 
conductance for both charge and energy transport deter- 
mine the difference and sum of Gec and Gca respec- 
tively. Thus the conductances of electron cotunncling 
and crossed Andreev reflection are experimentally acces- 
sible and can be compared to results from theoretical 
models. 



III. CIRCUIT THEORY 

The rules of circuit theory allows calculation of the 
cavity Green's functions in a network when the termi- 
nal Green's functions and structure of the connectors is 
determined. The terminals are characterized by known 
quasiclassical equilibrium matrix Green's functions Gi in 
Nambu-Kcldysh spaced The Green's functions depend 
on quasiparticle energy, and terminal temperature and 
chemical potential. The cavity Green's function in our 
model is denoted G c . The Green's functions are 4x4 
matrices including 2x2 Keldysh space and 2x2 Nambu 
space. For an explanation of our standard matrix no- 
tation, see Appendix [A] The general expression for the 
matrix current through the connector between terminal 



where 



4 + T, ( / ) ({G i ,G c }-2) 



(6) 



(7) 



pii' commutes with the Green's functions in Eq. © 
since it can be expanded in anticommutators of Gj and 
G c £l The matrix inversion in Eq. ([7]) can be performed 
analytically in Keldysh space due to the symmetries of 
these matrices. The spectral charge current can be ob- 
tained from the expression for the matrix current as 

2r,t = Tr|<73/| K ' ) | /8e, and the spectral energy current 

as = Tr|l| K) | /8e. The K superscript denotes the 
Keldysh matrix block of the current. 

Correlation between quasiparticles with opposite exci- 
tation energy from the Fermi surface is induced in the 
cavity due to Andreev scattering at the superconducting 
terminal. Cooper pairs transferred from the supercon- 
ductor into the cavity give rise to an electron with excita- 
tion energy E and a hole with excitation energy — E. The 
electron and hole quantum wave functions are initially in 
phase, but a relative phase will arise due to a small mis- 
match of the wave vectors* 1 ^ Their wave vectors are k = 
A:f\/1 ± E/E-p where hp is the Fermi momentum and Ep 
the Fermi energy. The relevant transport energy scale 
is the maximum of the temperature k^T and bias volt- 
age eV. The phase difference between the electron and 
the hole becomes Acp ~ 2Er/h w 2max(eV r , k^T^/h, 
where the dwell time in the cavity is r. The dwell time 
will be discussed in the next paragraph. We denote by 
Erh = H/(2t) the effective Thouless energy of the cav- 
ity. Let us consider the regime of negligible dephasing, 
characterized by max(eV, k^T) <C E^h- The presence of 
a superconducting terminal leads to prevailing electron- 
hole correlations since Acj> ~ ma.x(eV, fceT)/i?Th <C 1- 
In the regime of complete dephasing, on the other hand, 
max(eV, k^T) E^h, and initial many-particle phase 
correlation is lost since Acj) is finite and can only be 
described statistically. Thus the induced superconduct- 
ing correlations due to Andreev scattering are lost, and 
the wave functions of the quasiparticlcs in the cavity 
are not in phase. This dephasing effect is described in 
circuit theory by an additional terminal for "leakage of 
coherence"* 1 ^ Note that no charge or energy flows into 
this terminal. Circuit theory emerges from a discretiza- 
tion of the Usadel equation, and the dephasing term 
stems from the energy term in this differential equation. 
The expression for the matrix current due to dephasing 
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J D =ie 2 v V c E[a 3 ±,G c ] /h, 



(8) 



where v is the density of states in the normal state and 
V c is the volume of the cavity. 
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Let us now discuss the dwell time t defined above. The 
dwell time can be expressed as r = e 2 ^oKi?Totai- -^Totai 
is the total resistance to escape from the system, and it 
includes contributions from the contacts and diffusion. 
Diffusion is modeled by representing the diffusive region 
as a network of cavities connected by tunnel-like conduc- 
tors with resistance times area r = pd. Here p is the resis- 
tivity and d the lattice size in the discretized network?^ 
These connectors contribute to i?Totai • When diffusion is 
the dominating contribution to i?Totah the definition of 
the effective Thouless energy gives E Th = hD/(2L 2 ) in 
agreement with the continuum theory from which the cir- 
cuit theory is derived. D is the diffusion constant and L 
the typical length between contacts. E^ti is the relevant 
energy scale for the proximity effect in diffusive systems 
with negligible contact resistances^ In this paper, how- 
ever, we will consider the opposite limit that i?Totai is 
dominated by the contact resistances. Spatial variation 
of the Green's function inside the system is neglected, 
and we may discretize with only one cavity. The effec- 
tive Thouless energy is in this case £rh = fi/ {2e 2 v Q RL), 
where R is the sum of the interface resistances in paral- 
lel times area. The contact resistances induce an energy 
scale for dephasing similar to systems where diffusion is 
the dominant contribution to i?Totai- In Section lIVI we 
show in numerical calculations that the effective Thou- 
less energy is the relevant energy scale for the proximity 
effect in the cavity. 

We assume that inelastic processes in the cavity can be 
neglected since the characteristic time for inelastic inter- 
action is assumed to be much larger than the dwell time. 
The cavity Green's function is determined by demanding 
matrix current conservation at each energy. The sum of 
all matrix currents flowing into the cavity should vanish, 



the retarded "Usadcl equation" of the cavity: 



2e^ 
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=1,2,S 
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T^n^G- 

l J n 



i—voVcEasl, G c 
h 



0. (9) 



This equation determines the Green's function on the 
cavity, G c . The retarded and advanced components of 
G c can be parameterized in terms of one complex func- 
tion 6(E) as Gf = 03 cos(0) + a 1 sin(0) and G£ = 
— <7 3 cos(6>*) + &i sin(#*)i22 The definition of the Green's 
function implies that Re{cos(#)} is the normalized, en- 
ergy dependent density of states in the cavity, v(E)/vq. 

The Keldysh part of the Green's function is parameter- 
ized as G* = Gf^ h c — h c G^, where h c = lh^ :C + a^hT^. 
The normal terminals have Green's functions G\(2) = 
03T3 + (<73/il,i(2) + l^T,i(2))(n + 1^2), and the Green's 
function of the superconducting terminal is Gs = oil. 
We have assumed that any bias voltage eV <C A, where 
A is the gap of the superconducting terminal. Therefore, 
the only transport process into S is Andreev reflection 
since there are no accessible quasiparticle states in this 
terminal. 

The retarded part of matrix current conservation, Eq. 
(|9"|). gives an equation that determines the pairing angle 



.e 2 e 2 



rp{i) 



=5*2 + 



sin(6>) 



E 



- 1 n 



cos(0) = 0. 

*K ^2 + Ti S) (sin(0)-l) 

(10) 



The physical effect on the spectral properties from the 
various terms can be understood by comparing this equa- 
tion to the corresponding diffusion equation for a bulk 
superconductor . 33 ' 34 ! 35 This is given in Eq. (II. 29b) of 
Ref. and becomes in our notation^ 

HD d 2 6 ( h \ H 

— t— 7+ \E - - — sinfl + Acos^ sin0cos6» = 0, 

2 dx 2 V 2t e J 

where D is the diffusion constant, A the gap to be deter- 
mined self-consistently, 1/te the inelastic scattering rate, 
and l/r s f the spin-flip scattering rate. Comparing this to 
Eq. (flQ|) , we see that the coupling to the superconductor 
induces superconducting correlations, and that the cou- 
pling to the normal terminals give quasiparticles a finite 
lifetime. Spin-flip scattering could be included by tak- 
ing into account magnetic impurities in the cavity. We 
consider a normal metal cavity. To describe a supercon- 
ducting cavity, we would have to include a pairing term 
in the Hamiltonian. This would result in a term in Eq. 
([9]) with the same structure as the term corresponding to 
coupling to the superconducting terminal. As long as we 
do not consider Josephson effects, the effect of supercon- 
ductivity in the cavity could therefore be included by a 
quantitative renormalization of the coupling strength to 
S, which is straightforward. 

In the regime of complete dephasing, the term propor- 
tional to sin(0) dominates in Eq. (TIT))) because of the 
large factor E/E^. The solution in this limit is = 0, 
which means that there are no electron-hole correlations 
in the cavity. 

Let us now consider the Keldysh part of Eq. ([9]). We 
take the trace of this matrix block after first multiply- 
ing it with 03. The resulting equation determines the 
distribution function /it,c hi the cavity, 

Gt,1 (h>T,l — ^T,c) + Gt,2 (h>T,2 ~ ^T,c) 

+ G T ,s (0 - hr, c ) = 0. (11) 

This implies charge conservation at each energy, with ef- 
fective energy dependent conductances Gt.i between the 
cavity and terminals Ni, N2 and S. The zero in the term 
Gt,s (0 — ^t,c) represents the charge distribution func- 
tion in the superconductor, /it,S; which vanishes since 
the superconductor is grounded. The conductance Gt,s 
controls the Andreev reflection rate in the cavity since it 
is the conductance between the cavity and the supercon- 
ducting terminal. The effective conductances are given 
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in terms of the pairing angle and the transmission prob- 
abilities as 



(2 - T^)Re {cos(0)} + T^D T 



nn *— ' 



2 + T^ (cos(0) - 1) 



(12a) 



2 + Ti S) (sin(0)-l) 



(12b) 



Here L» T = [Re {cos(<9)}] 2 + [Re {sin(6»)}] 2 and i = 1,2. 
The term proportional to Gt,s m Eq. (jlip describes 
conversion of quasiparticlcs in the cavity into conden- 
sate in the superconducting terminal. There is an anal- 
ogous term in the Boltzmann equation for a continuum 
superconductor; 35 ' 36 which describes conversion between 
quasiparticles and superconducting condensate over the 
coherence length. This analogous term has a similar de- 
pendence on 9. The rate of this conversion is controlled 
by A in the continuum case, and by Gt,s m the cavity 
of our discretized theory 

In the regime of complete dephasing 8 = 0, 
the conductances to the normal terminals Gx,i coin- 
cide with the Landaucr-Buttiker formula. Gt,s = 



2e 2 J2 n (Tk a ') 2 /(TrH(2 - T.P) 2 ) corresponds to the An- 
dreev conductance of an N/S interface in a diffusive sys- 
tem, calculated by Beenakker.— Thus the distribution 
function /it, c can be determined in the regime of com- 
plete dephasing from well-known results by demanding 
charge current conservation. 

In the tunnel barrier limit, all transmission proba- 
bilities are small, and we can expand to first order in 



Tn in Eqns. (p]). We define g t = e 2 Tk l '/(irH) for 
i = 1,2, S. For the normal terminals (i = 1,2) we find 
Gt,? = 5?Re {cos(#)}. Re{cos(6>)} gives the normalized 
density of states in the cavity which is under the influ- 
ence of the proximity effect. The tunnel conductance to 
the superconductor becomes Gt,s = SsRe {sin((9)} which 
vanishes when there is complete dephasing (9 = 0. This 
is expected since the Andrecv conductance of a tunnel 
barrier between incoherent normal and superconducting 
terminals vanishes^ 

The trace of the Keldysh block of Eq. ((9]) gives an 
equation that determines the distribution function /il, c 
in the cavity, 



G L ,1 (h L>1 - /l L ,c) + G L , 2 (/l L ,2 - h L<c ) = 0. 



(13) 



This is energy conservation at each energy, with effec- 
tive energy dependent conductances for energy transport 
Gl,;- No energy current can flow through the contact be- 
tween the cavity and the superconducting terminal since 
no net energy is transferred into S by Andreev reflection. 
Our calculation is restricted to E <C A, but in the general 
case a quasiparticle current which carries energy can flow 



into the superconducting terminal for E > A. The effec- 
tive conductances for energy transport are given in terms 
of the pairing angle and the transmission probabilities as 



(2-Ti i) )Re{cos(9)}+T^D 



W; 



2 + (cos(0) - 1) 



(14) 



Here D L = [Re {cos(6>)}] - [Im {sin(0)}f and i= 1,2. 

In the tunnel barrier limit we find that Gl,z = Gx,j = 
t/iRe {cos(0)} for £ = 1,2, which means that the effective 
conductance for energy transport and charge transport 
into the normal terminals are equal. The conductances 
correspond to the usual quasiparticle tunnel conductance 
in this case. 

Eqns. (TIT))) . (TIT]) and (TiUj) determine the Green's func- 
tion in the cavity. The charge and energy currents, Jr,i 
and between terminal i and the cavity can be cal- 
culated once G c is known. Comparison of Eqns. j2| and 
(|4"| with the expressions for 7t.i and 7l.i obtained from 
circuit theory, allows us to determine the conductances 
associated with the various transport processes: 



G a(E) — - [ < ~* T ' 1 ^ T ' 2 — Gt - s ) Gl,iGl,2 
4 \Gt,i + Gt,2 + Gt,s Gl,i + Gl,2 



(15a) 



„ ,„> 1 / Gl,iGl,2 Gt,iGt,2 

Gec(^) = 9 7; — Tr — ± 



Gt i + G 



T,2 



G 



T.S 



(15b) 



These formulas are the main result of our calculation. 
Eq. (|15b[) shows that Gec — Gca is positive. Thus the 
nonlocal conductance dI\ldVi is positive. In the limit 
where the coupling to the superconducting terminal van- 

ishes, i.e. all T„ — > 0, only Gec remains nonzero and 
the conductance agrees with the result for a normal dou- 
ble barrier system. If the conductance to one of the nor- 

(i) 

mal terminals vanishes, i.e. all T„ — > for e.g. i = 2, 
only Gda is nonzero. When the coupling between the 
superconducting terminal and the cavity is very strong 
Gt,s S> Gt.i, Gt,2, we recover the result of Ref. [H that 
Gec = Gca, which means that the nonlocal conduc- 
tance vanishes since transport by electron cotunncling is 
exactly canceled by crossed Andreev reflection. 

To describe a device where spatial variation in a bulk 
region is important, a model with several cavities con- 
nected in a network is required. The connectors between 
cavities represent the intrinsic resistance due to diffusion, 
and contribute to i?Totai and thus to the effective Thou- 
less energy E^h = VO^^oK-RTotai)- However, as long 
as there are no Josephson currents in the network, the 
symmetry Gec(ca) = a + ( — )P of Eqn. (|15b[) persists. 
Here a and /3 are positive numbers. This follows since 
the currents flowing out of the normal terminals given by 



circuit theory can be written I^^E) 
and It^(E) 



Cti%\i 



G[ l) ( Vl - &L.2) 

C™h r r,2 where Cj/ , and 
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(7^2 are coefficients. Comparing these expressions to 
Eqns. ^ and dU), we see e.g. that Gec — Gca = 
2eC^l > for contact i, regardless of the internal struc- 
ture of the network of cavities. Thus the sign of the non- 
local conductance is not affected by diffusion or network 
geometry. We interpret this result as the consequence of 
a symmetry between crossed Andreev reflection and elec- 
tron cotunncling. Both processes involve the transfer of 
quasiparticles through the contacts to the normal metals 
and the network between them, but crossed Andreev re- 
flection also involves Andreev reflection at the interface 
to the superconducting terminal. Thus the resistance 
limiting crossed Andreev reflection can at minimum be 
as small as the resistance for electron cotunncling unless 
other physical processes affect these quantities. We be- 
lieve that an explanation of the measurements in Ref. 
requires additional physical effects not considered here. 



A. Analytically solvable limits 

In this section, we give results for two limits where it is 
possible to solve Eq. (TIT))) analytically and obtain simple 
expressions for the conductances in Eqns. (|15[) . 



1. The regime of complete dephasing 

When E 3> £rh there is complete dephasing. In this 
case the conductances Gt,, = Gh,i = g% for i = 1,2 
agree with the Landauer-Biittiker formula and Gt,s is 
Bcenakkcr's result for the conductance of an N — S inter- 
face in a diffusive system, as noted above. In this limit, 
the conductances in Eq. ([15)) become 



G 



G 



DA 



ec(E) 

CA 



9l 



Gt,: 



4 (ffi + 32) (.9i + .92 + G T ,s) 

9\9i 1 

2 (gi +52H51 +92 + G T ,s) 

< f2(. gi +.g 2 ) + GT,s for EC 
* ] G T ,s for CA 



(16a) 



(16b) 



The currents resulting from these expressions using Eq. 
(fl]) can also be calculated from a semiclassical approach 
by demanding charge conservation in the cavity using the 
well-known theory for incoherent N/S transport. The re- 
sult for Gec — Gca was shown in Ref. 0, where also 
spin polarizing contacts are employed to give a negative 
nonlocal conductance. Again, in the limit of strong cou- 
pling to the superconductor, ^/Gt.s < 1, we see that 
Gec = Gca and thus the nonlocal conductance vanishes. 



2. Tunnel barrier limit 



In the tunnel barrier limit, all transmission probabil- 

(i) 

ities are small. We expand to first order in T„ in the 



expressions for the matrix currents. This corresponds to 
putting p$ -f 1/4 in Eq. ©. The resulting Eq. © can 
be solved analytically without resorting to a parametriza- 
tion of G C J£ The solution is of course equivalent to what 
we obtain from Eqns. (jTOJ) , (fTTj) and (fl3|) in the same 
limit. Let us consider the resulting expressions for the 
conductances of nonlocal transport in some limits. At 
zero energy, we obtain 



Gda — 
Gec = 
Gca = 



g\ 



gi 



4 [(3i+52) 2 +5l] 3/2 ' 

g\gi 2(gi+g 2 ) 2 + g| 

2 [(.9i+.9 2 ) 2 + g§] 3 / 2: 
5152 



95 



2 [{gi+g2? + glYl 2 ' 



(17a) 
(17b) 
(17c) 



These results correspond to completely phase coherent 
transport. The full counting statistics for the same sys- 
tem in this regime has been calculated in Ref. HH. In 
that paper, it is found that the cross correlation noise 
can have both signs in the three terminal device. 

If there is complete dephasing, the cavity spectral 
properties are like those of a normal metal. This gives 
Gt.s = since the Andreev conductance of a tunnel bar- 
rier vanishes for incoherent N/S transport^ Therefore, 
there are no transport channels into the superconduct- 
ing terminal and it is effectively isolated from the cir- 
cuit. Only the conductance for electron cotunneling is 
nonzero, and we obtain from Eq. (|17b[) 



Gec 



1 



1/.91 + V.92 



(18) 



i.e. addition of the conductances between the cavity and 
both normal metal terminals in series. This result corre- 
sponds to normal state tunneling between Ni and N 2 . 



IV. RESULTS 

In an experimental situation, nonlocal transport can 
be probed in measurements of the voltage or current 
in terminal Ni resulting from the injection of current 
through terminal N 2 i^ At zero temperature, the nonlo- 
cal differential conductance 9-f c harge,i/^^2 as a function 
of eV 2 / i?Th corresponds to Gec — Gca , see Eq. ([3]) . This 
quantity is given by Eqns. (|15bp and is always positive 
within our model. The detailed behavior as a function of 
eV 2 depends on the pairing angle of the cavity, 9, which 
is determined by Eq. (flO|) . This equation can be solved 
numerically, and determines the spectral properties and 
thus the nonlocal conductances when the parameters of 
the model are fixed. We will show results for three sys- 
tems with different combinations of contacts which can 
be fabricated with state of the art nanotcchnology. The 
contacts we study are tunnel barriers or metallic con- 
tacts. The systems represent various realizations where 
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our circuit theory model applies. We show that the non- 
local conductances are very sensitive to the type of con- 
tacts in the system, and have a strong dependence on 
the Thouless energy. We have also investigated the effect 
of inelastic scattering inside the cavity, and have found 
no notable qualitative differences on the conductances as 
compared to the elastic case, thus these results are not 
shown here. 



A. Tunnel barriers 

When all the three contacts are tunnel barriers, the 
equation for matrix current conservation can be solved 
analytically as mentioned in Section riff A 21 The result- 
ing expressions are quite complicated. We expand to first 

(i) 

order in the transmission probabilities since Tk > < 1 for 
i = 1,2, S, and let <?i denote tunnel barrier conductance 
in the normal metal state. Let us first consider a symmet- 
ric system, gi/gs = 92/ gs = 0.1. We define the Thouless 
energy in this case as E^h = %s/(2e 2 ^oT4). The spec- 
tral properties of the cavity are shown in Fig. [3J In 
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FIG. 3: (Color online) Spectral properties of the cavity when 
all contacts are of tunnel type with gi/gs — 32/gs = 0.1. 
Solid line (blue) shows normalized density of states, u(E)/vo- 
Dashed line (red) shows Gt,s /gs which is the parameter that 
determines the Andreev reflection rate of quasiparticles. 

this plot we show the normalized density of states in the 
cavity and the conductance Gt,s which controls the An- 
dreev reflection rate. From Fig. [3] we see that below 
the Thouless energy of the cavity, the quasiparticle den- 
sity of states is suppressed due to electron-hole correla- 
tions. This affects all transport processes since they rely 
on quasiparticles propagating through the cavity. Above 
the Thouless energy, the density of states approaches the 
value in the normal state. This is the typical behavior for 
proximity coupled systems^ The conductance between 
the cavity and the superconducting terminal, Gt,s, ap- 
proaches the normal state value at low energies. At high 
energies, Gt,s vanishes, as expected for tunnel barriers 
in the regime E < A when the proximity effect can be 



neglected This affects crossed and direct Andreev re- 
flection, which vanish when Gt,s goes to zero. 

A plot of the conductances for electron cotunneling, 
and crossed and direct Andreev reflection is shown in 
Fig. 3J The conductances have a rapid increase near 
E = Et\i related to the energy dependence of the den- 
sity of states. At high energies Gec approaches the value 
for normal double barrier tunneling, Eq. ((THJ) . Gqa and 
Gda vanish at high energy. These conductances are de- 
termined by an interplay of the density of states and 
the Andreev reflection rate, and are small when cither 
of these quantities are small. The measurable nonlocal 
differential conductance, Gec — Gca, is a monotonously 
increasing function with increasing energy, which starts 
at small value and approaches Gec above the Thouless 
energy. 

0.06 1 1 1 1 . 1 




FIG. 4: (Color online) Conductance for the various transport 
processes through the cavity when all contacts are of tunnel 
type. The two normal contacts are have the same conductance 
in this plot, gi/gs = 52/gs = 0.1. Solid line (red) Gec/ss, 
dashed line (blue) Gca/<?s, dash-dot line (green) Gda/^s, and 
dotted line (purple) (Gec — Gca)/ps- 

Let us now consider the effect of asymmetry between 
the tunnel barriers to the normal metals. The expres- 
sions for the conductances at zero energy, Eqns. fll7[) . 
show that the direct Andreev conductance of one con- 
tact, Gpj^ for example, is proportional to g\ since two 
electrons have to tunnel through the connector with con- 
ductance g\. G^ A is, however, only weakly dependent 
on g2 ■ The same is true for the direct Andreev conduc- 
tance of connector 2 with 31 <-» 32- The direct Andreev 
conductances are therefore relatively independent on the 
asymmetry On the other hand, nonlocal conductances 
are proportional to 31(72 because a quasiparticle has to 
tunnel through both connectors. These conductances are 
sensitive to the asymmetry which we define as a = 31/32 ■ 
For the conductance measured at terminal Ni we see that 
^ec(CA)/^da 01 l/ a - Conversely, the conductance mea- 
sured at terminal N2 gives Gj 3( i ( - CA j/Gj- 1 ^ oc a. Thus 
asymmetry suppresses the nonlocal conductance of one 
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contact, and enhances the nonlocal conductance of the 
other contact. In Fig. [5] we show the conductances for a 
system where gi/gs = 0.1 and = 0.3. The asym- 

metry is now a = 1/3, thus the effect of nonlocal pro- 
cesses is enhanced when we consider the conductances 
measured at terminal Ni. The spectral properties in this 
case are similar to those of the symmetric system shown 
in Fig. [31 and are not shown here. Comparing Figs. [31 
andfSJwe see that the conductances for crossed and direct 
Andreev reflection are not as peaked in the asymmetric 
system as in the symmetric system. We also see that 
the relative magnitude of the nonlocal conductances to 
the direct Andreev conductance has increased. In the re- 
mainder of the manuscript, we consider the conductances 
of contact 1 in asymmetric structures a = 1/3 since we 
are mostly interested in the conductance resulting from 
nonlocal processes. 
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FIG. 5: (Color online) Conductance for the various transport 
processes through the cavity when all contacts are of tunnel 
type. Conductances are normalized by gs, and parameter 
values are gi/gs = 0.1, 32/gs = 0.3. Solid line (red) Gec/<?s, 
dashed line (blue) Gca/3s, dash-dot line (green) Gda/<?s, and 
dotted line (purple) (Gec - Gca)/5s- 



B. Tunnel contacts to Ni and N2, 
metallic contact to S 

Systems where two normal metal terminals are con- 
nected to a cavity, in which the cavity may itself be part 
of a larger superconducting structure, can be studied in 
our model. The contact to the superconductor could in 
this case be e.g. diffusive or metallic. In this section, 
we assume tunnel contacts to the two normal metals, 
and metallic contact to the superconductor. The metal- 
lic contact is described by the transmission probabilities 

(s) 

Tn = 1 for all propagating channels and zero otherwise. 
We choose parameters gi/gs = 0.1 and 32 /.9s = 0.3. The 
Thouless energy is £rh = %s/(2e 2 i/oVQ. In Fig. [5] we 
show the spectral properties of the cavity. Below the 
Thouless energy, the quasiparticle density of states is sup- 



pressed. At low energy, Gt,s approaches the value in the 
normal state. This is similar to the reentrant behavior 
in diffusive systems i 31 i 39 With increasing energy, the rate 
of Andreev reflection at the superconducting terminal in- 
creases. The value of Gt,s/<?s reaches its maximum value 
2 for E/Et\i ~ 20. This result agrees with the Blonder- 
Tinkham-Klapwijk formula for high transmission proba- 
bilities, valid when proximity effect is negligible. 38 




FIG. 6: (Color online) Spectral properties of the cavity when 
the contact to the superconducting reservoir is metallic with 
TA ' = 1 for all propagating channels and zero otherwise, and 
the contacts to the normal reservoirs are of tunnel type with 
5i /.9s =0.1 and gi/gs = 0.3. Solid line (blue) shows normal- 
ized density of states in the cavity, v(E)/vo, and dashed line 
(red) shows Gt,s/<?s which is the parameter that determines 
the Andreev reflection rate of quasiparticles in the cavity. 

The conductances for this system are shown in Fig. [7] 
All conductances are suppressed below the Thouless en- 
ergy by the low density of states. Above the Thouless 
energy, the increasing Andreev reflection rate leads to a 
suppression of Gec an d an enhancement of Gqa- The 
overall behavior of the nonlocal conductances is deter- 
mined by the interplay of the dependence on the density 
of states and the Andreev reflection rate. The two con- 
ductances associated with Andreev reflection, Gca and 
GdA; reach their maximum at approximately E / '£rh = 4 
and then decrease slowly, reaching their final value at 
E/E Th w 20. 

C. Metallic contacts to Ni and N2, 
tunnel contact to S 

Let us now consider a system where the normal ter- 
minals Ni and N2 are connected to the cavity through 
metallic contacts described by T„ = 1 for all propagat- 
ing channels and zero otherwise. The superconducting 
terminal is connected by a tunnel barrier of conductance 
gs, and the Thouless energy is £rh = %s/(2e 2 ^oV c ). 
The spectral properties of this system are very similar to 
the tunneling case in Section llV A[ and we do not show 
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FIG. 7: (Color online) Conductances for the various transport 
processes when the contact to the superconducting reservoir is 

(S) 

metallic with T„ — 1 for all propagating channels and zero 
otherwise, and the contacts to the normal reservoirs are of 
tunnel type with gi/gs =0.1 and gi/gs — 0.3. Solid line (red) 
Gbc/ss, dashed line (blue) Gca/<?s, dash-dot line (green) 
GoA/gs, and dotted line (purple) (Gec - Gca)/3s- 



them here. However, the conductances between the nor- 
mal terminals and the cavity are qualitatively different 
in the present case. In Fig. [5] we show Gt.i for both the 
tunneling case in Scction [lV Al and where the normal ter- 
minals are connected by metallic contacts. At energies 
below the Thouless energy, Gt,i is in the tunneling case 
qualitatively similar to the density of states. However, 
with metallic contacts to the normal reservoirs, Gt,i is 
large at zero energy and decreases as the energy increases 
beyond £rh- At high energy, Gt,i is equal for the two 
cases and corresponds to the result in the normal state 
since gi/gs is the same for the two curves. Gt,i for a 
metallic contact is qualitatively similar to the conduc- 
tance of a metallic normal metal- superconductor inter- 
face, except that in this case the "superconductor" is the 
cavity which is under the influence of the proximity effect 
from S and the relevant energy scale is i/ph instead of A. 

The conductances for the present system are shown in 
Fig. [5J A new feature here is a small dip in Gec at 
E = i?Th before the rapid increase above Erh- Gec is 
proportional to Gt,i and inversely proportional to Gt,s- 
The dip in Gec can therefore be understood by the de- 
creasing charge current conductance (Gt,i) and the peak 
in Gt.s around E = i?Th- At higher energy Gec in- 
creases as Gt,s vanishes. The dip in Gec leads to a 
larger dip in Gec — Gca since Gca increases with in- 
creasing energy at E < i?Th- In the tunneling case of 
Section IIV Al the increase of the Andreev conductance 
is compensated by an increasing charge current conduc- 
tance in Gec — Gca, and there is no dip in the nonlocal 
differential conductance. At high energy, Gca and Gda 
vanish since the Andreev reflection rate vanishes. 



FIG. 8: (Color online) Energy-dependence of conductance 
Gt,i for the present case of metallic contacts to normal reser- 
voirs and tunnel contact to the superconducting reservoir in 
solid line (blue), and for the case with only tunnel barriers 
studied in Sec. IIV Al in dashed line (red) . In both cases we 
have put the parameters gi/gs = 0.1, and g?/gs = 0.3 



0.08 




FIG. 9: (Color online) Conductances for the various trans- 
port processes when the contacts to the normal reservoirs are 
of metallic type with = 1 where i = 1,2 for all propa- 
gating channels and zero otherwise, and the contact to the 
superconducting reservoir is of tunnel type with conductance 
gs- In this plot gi/gs — 0.1 and <?2/<7s = 0.3. Solid line (red) 
Gec/<?s, dashed line (blue) Gca/<7s, dash-dot line (green) 
Gda/ps, and dotted line (purple) (Gec - Gca)/.9s- 



V. CONCLUSION 

In conclusion, wc have studied nonlocal transport in a 
three terminal device with two normal metal terminals 
and one superconducting terminal. To this end we have 
applied the circuit theory of mesoscopic transport. The 
connectors between the circuit elements are represented 
by general expressions, relevant for a wide range of con- 
tacts. Dephasing is taken into account, and is governed 
by the inverse dwell time. This gives rise to an effective 
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Thoulcss energy. For this model we have calculated the 
conductance associated with crossed Andreev reflection 
and electron cotunneling, as well as direct Andreev reflec- 
tion. We find that for any contacts, electron cotunneling 
is the dominant nonlocal transport process. Results for 
various combinations of contacts of tunnel and metallic 
type are shown, and this demonstrates the strong depen- 
dence on the nature of the contacts. The conductance 
of nonlocal transport has a strong dependence on the 
Thouless energy. This is because at these energies, de- 
phasing leads to a loss of induced superconducting corre- 
lations. Several of the characteristics of our model have 
been observed in experiments, although we do not make 
complete qualitative contact to experimental data since 
electron cotunneling always dominates crossed Andreev 
reflection. We suggest determining the conductance due 
to these processes independently by carrying out addi- 
tional energy transport measurements. 
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APPENDIX A: NOTATION 

Matrices in Nambu and Keldysh matrix space are de- 
noted by a "hat" (M) and a "bar" (M) respectively. The 
symbol used for a unit matrix is 1, and the Pauli matri- 
ces are denoted a n and f„ in Nambu and Keldysh space 
where n = 1,2,3. Compositions of matrices in Nambu 
and Keldysh space are formed by a direct product to 
make up 4x4 matrices in Nambu-Kcldysh matrix space, 
so that e.g. 



03 n = 



'0010 
-1 
10 

vO -1 



(Al) 
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The distribution functions /it(-E), hi,(E) can be written in 
terms of the particle distribution function f(E) as fox = 
1 - f(E) - f(-E) and h L = -f(E) + f(-E), see Ref. M- 
To make contact with our notation we must identify a = 
cos and j3 = sin# in Eq. (II. 29b) of Ref. HI In addition to 
this, we put 7 = since there is no such component in the 
retarded Green's function of the cavity, <E> = since we are 
considering stationary phenomena and vs = since there 
is no supercurrent in the cavity. 



